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ACTIVITIES (January - June 1993): 

1. Experiments with thermal imaging infrared radiometers. A new AGEMA 900 series 

radiometer, equipped with an array of filters and calibrated, was purchased by Michigan 
Tech. Our group did field experiments and designed a laboratory system for controlled 
study of simulated ash clouds in our laboratory. ^ 

2. Participation in the Federal Aviation Administration Workshop on Old Volcanic Ash 
Clouds in Washington, D. C. Presented two papers (see attached). 

3 Development of radiative transfer method to retrieve particle sizes, optical depth, 
and particle mass burdens in volcanic clouds, using AVHRR thermal infrared bands 4 
and 5. Finished and submitted paper to JGR/Atmospheres, completing one of the 
proposals principal objectives (see attached). 


ACTIVITIES (July - December 1993): 

1. Development of a radiative transfer model to retrieve particle sizes, optical depth and 
volcanic ash mass burdens in volcanic clouds, using TIMS, MODIS and ASTER. The 
new generation of sensors will allow us to apply the same principles to multiple bands 
of data, allowing better constrained and more robust results. 

2. Participation in American Geophysical Union Fall Meeting. Presentation of three 
papers. 
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demonstrate spectral changes associated with morning solar 
warming, showing that there are large changes in spectra over 
the period of 8-10 AM local time when TM data are typically 
collected over Central America. 
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presentation will illustrate the problems of automation wi 
examples of AVHRR volcano images studied in detail. 
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Satellite-Baaed Estimates of Ash Content of Drifting Volcanic Clouds 

W I Rose (G eo logical Engineering, Michigan Tech Univ, Houghton, MI 49931; 906 487-2531; 
e-mail: iaman@mtn.edn) S. Wen, D. J. Schneider 

Over the past year we have developed a radiative transfer model that can retreive particle sizes, 
optical depth and particle masses in transparent volcanic clouds (Wen and Rose, this meeting) 
iigitig two of infr ared data from the AVHRR. The method is a new tool to 

wamiiw the fate and transport of volcanic ash in earth’s atmosphere. The two band AVHRR 
method can be applied to a 15 year archive of data on volcanic clouds beginning in 1978. We 
expect that the results will help us measure the total volumes of eruptions, by facilitating 
measurement of far flung ash. The study of multiple images from the same eruption will enable 
study of dispersion and particle fallout over the first days of the drifting cloud (Schneider et al, 
this meeting). More robust methods expanding the same ideas to multiple thermal infrared 
bands are being developed using the TIMS, MODES and ASTER sensors. These methods will 
be overdetermined and will allow us to refine the techniques considerably. We have also begun 
study of simultaneous data on volcanic clouds collected by AVHRR and TOMS, using processed 
data obtained from the TOMS investigators at NASA Goddard (Hossli et al, this meeting). 
Because the AVHRR detects ash, and the TOMS detects SO* using both sensors potentially 
allows us to see the separation of ash and SO* 
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of the cloud. We hope to invert the bispectral AVHRRdata to estimate 
concentration. 
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RETRIEVAL OF PARTICLE SIZES AND MASSES IN 
VOLCANIC CLOUDS USING AVHRR BANDS 4 AND 5 


Shiming Wen and William L Rose 


Department of Geological Engineering, Geology and Geophysics 
Michigan Technological University 
Houghton, MI 49931 


Submitted for publication to JGR- Atmospheres: August 23, 1993 



Retrieval of Particle Sizes and Masses in 
Volcanic Clouds Using AVHRR Bands 4 and 5 


Shiming Wen and William L Rose 

Department of Geological Engineering. Geology and Geophysics 
Michigan Technological University 
Houghton, MI 49931 


ABSTRACT 

The AVHRR sensor on polar orbiting NOAA satellites can discriminate between 
volcanic clouds and meteorological ones using two-band data in the thermal infrared. 
•Hus paper is aimed at developing a retrieval of the particle sizes, optical depth and 
particle masses from AVHRR two-band data of volcanic clouds. Radiative transfer 
calculations are used with a nansparent cloud model that is based on assumptions of 
spherical particle shape, a homogeneous underlying surface and a simple thin cloud 
parallel to tire surface. The model is applied to observed AVHRR dam from a 13 hour old 
drifting cloud from tire 19 August 1992 eruption of Crater Peak/Spmx Volcano. Alaska. 
The AVHRR dam fit in the range of results calculated by the model, which supports its 
credibility. According to the model results, the average of effective particle radius in the 
test frame of tins cloud is in the range of 2 to 2.5 pm, the effective emissivity averages 
0 6 and the optical depth is about 0.60 - 0.65. Hie mass of ash estimated amounts to 
35 000 -50.000 tons in the test frame of the cloud, and to 0.20 - 0.28 x V? tons in the 
whole cloud. The total estimated mass is about 0.5-0.891 of the mass measured m die 
blanket Applications of the new retrieval method are listed. 


INTRODUCTION 


is 


Measuring the size and burden of silicates and other components in volcanic clouds 
of interest to those studying volcano atmosphere interactions (Table 1). Weather 



.■...lli.-i have been a useful way to track drifting volcanic clouds (Hanstrum and Watson, 
1983; Sawada, 1987) and two-band processing of thermal infrared data from weather 
has allowed for the discrimination of volcanic and meteorological clouds 
(Prata,1989a; Holasek and Rose, 1991). The explanation of the cause of infrared two-band 
and volcanic clouds have been discussed by Prata (1989a, 1989b), Holasek 
and Rose (1991) and Schneider and Rose (1993), and is the result of scattering and 
absorption of thermal emission from matter underneath the volcanic cloud by the cloud 
itself. In this paper we develop a model of radiative transfer to attempt to retrieve the 
particle sizes and masses of particles in drifting volcanic clouds. The model builds on 
work of Prata (1989b) on volcanic clouds, on methodology published by Yamanoucht et 
al (1987), and on cloud retrieval methods of Lin and Coakley (1993). We compare our 
method with actual data on the Crater Peak/Spurr eruption of 19 August 1992 to begin to 
evaluate the model, and discuss some of the uncertainties and applications. 

BASIC THEORY 

Radiative transfer calculations have been used to develop two-band models for 
retrieving the optical parameters of clouds, such as particle sizes, emissmty, 
transmissivity, and cloud cover, because radiance attenuation through the atmosphere 
functionally depends on geometrical and optical properties of the clouds in the process of 
radiative transfer. The observed radiances by a satellite-based remote sensor through a 
transparent cloud is composed of two parts, i.e., radiance from the clouds and from the 
underlying surface. Generally, if the fraction of partial cloud cover in a field of view is 
taken into account (Coakley, 1983; Lin and Coakley; 1993), a linear model is valid under 
the following assumptions; (1) the cloud approximates a planar homogeneous cloud layer 
parallel to the surface (or a single-layer cloud system); (2) the background surface is 
homogenous; and (3) the atmosphere above the cloud and between the surface and the 
cloud are dear windows. In this case the observed radiance /,in a narrow band < centered 
at wavelength \ is given by the following equation: 



J 1 -(1-A C )B(T.) +A e (2’iB(T c ) -*-t i B(r,) ) 


( 1 ) 


where T, is the brightness temperatures of the surface, T c the temperature of the top of 
the cloud, B the Plank function, A c the fraction of the clouds in the field of view, e, the 
emissivity, and transmissivity of the clouds. The pixels are partially covered by clouds 
if a c is less than 1. If the clouds are optically thick and completely overcast (constant 
and A= 1), the measured radiance /, approximately equals to B(TJ. On the other hand, for 
a partially transparent cloud layer overlying a warm surface ( B(T t ) » B(TJ\ the radiance 
(I) approaches as % a measure of the transparency, approaches 0. 


Radiative Transfer Equation 


To obtain the theoretical radiance defined in equation (1), the radiative transfer 
calculation has to be used because of unknown emissivity and transmissivity at different 
wavelengths. The radiative transfer equation for a scattering plane parallel atmosphere is 

given by: 

u.^£ <1>) (2) 

oX 


J(T,p,<j>) is the source function given by 




T5tJo J-i 


where I is the diffuse radiance, 7 tF G incident solar flux, T the optical depth, p the cosine 
of the zenith angle, co the single scattering albedo is defined as the radio of scattering 
cross section a, to extinction (scattering plus absorption) cross section a 6 , p 0 the cosme of 
the solar zenith angle, and PQi,< fai’.f) is the axially-symmetric phase function definmg 
the light incident at p\ f which is scattered in the direction p, ty. The smgle-scattering 
albedo, phase function, and extinction cross section are functions of the incident 
wavelength, particle size and shape, and refractive index. In the 10-12 pm region, the 


source due to incident sunlight, i.e. the last term S^IWoFqPOi^-.Po^o) 6 '^ subscript 0, 
is negligible compared with radiation due to emission. 


Eddington’s Approximation 


To solve equation (2), several methods have been developed, such as Discrete- 
ordinates method (Chandrasekhar, 1960; Liou, 1973), Adding method (van de Hulst, 
1980), Monte Carlo method. Doubling method. Two-stream and Eddington s 
approximation (Eddington, 1916). Here Eddington’s approximation was used because it 
is a good approximations for thick layers (Liou, 1992) and offers a way to rapidly 
compute irradiances with an accuracy of several percent (Shettle and Weinman, 1970). 
Eddington’s approximation assumes that the radiance can be simply approximated by a 

linear function 

J(T,p)-r < ,(T)+J x (‘e)|i (3) 

and the phase function can be approximated by 

where the asymmetric parameter g=(l/2)JP(0)cos0d(cos0) (0 is the angle between incident 
and scattered radiances) is a measure of the amount of radiation scattered toward the 

Earth’s surface relative to that scattered back to space. 

If an atmosphere composed of homogeneous layers is assumed, we have 
3(o/3x=3g/9T=0, and straightforward analyses yield the following solution: 

j »Ae' fcT +Be* kT (5) 

O 



3 (l-co) 3(i-a» Be .K 

l-cog N T^cog 


( 6 ) 


Therefore, the downward (F ) and upward flux are given by 



F-(T)»2£nJ(^)d^=-(l+C7)Ae- ,rt -(l-C7)Be*‘ T (7) 

FMX)»2 < £nJ(^)djl*(l-C7)Ae* J "+(l+C7)Be^ (8) 

where Jfc=[(l-co)(l-oog)] w and C/=[4(l-co)/3(l-cog)] . 

For the boundary condition for a homogenous parallel cloud, it is assumed that the 
medium is only illuminated from above by a known source of radiation: 

F" ( 0 ) S ~F 0 , F* ( T c ) *0 . 

From this boundary condition, coefficients A and B can be determined. 

From above, the reflectivity of the atmosphere (R) is the proportion of reflective flux 
at the top of the atmosphere to the incident flux, and the transmissivity of the atmosphere 
is the proportion of transmitted flux to incident flux, i.e., 

- F*{0) _ (1-L7 2 ) (e^-e**-) 

F‘(T C ) 4 U (10) 

' T a (l-U) 2 e _Jrt — (l+U) 2 e’ J ”* 

and emissivity E=l-R-T. 

TEMPERATURE DIFFERENCE MODEL 

Prata (1989b), Wu (1987), and Yamanouchi et aL (1987) applied the temperature 
difference method to retrieve weather cloud parameters for complete cloud cover. In fact, 
the radiative calculation for the temperature difference method is a special case of 
equation (1) under the assumption A=l. The study frame selected (see below, Fig. 4) is 
located in the center of cloud image, and the pixels of the frame can be reasonably 


considered as complete cloud cover, but still partially transparent. 

Band 4 (10.3 -11.3 pm) and band 5 (11.5 - 12.5 pm) of AVHRR are used in our 
temperature difference model, because the temperature difference of band 4 and band 5 
(T 4 - Tj) is a fairly reliable criterion used to distinguish meteorological clouds and volcanic 
clouds. When a volcanic cloud exists, T 4 - T s is usually < 0, otherwise T 4 - T s > 0 (Prata, 
1989a; Schneider et al., 1993). 

The purpose of the Temperature Difference Model (TDM) proposed here is to 
estimate the mass of particles in volcanic clouds. So two parameters, particle radius and 
optical depth, are essential in this model. The radiative transfer model is based on the 
following assumptions: (1) the shape of the volcanic particles is spherical, therefore Mie 
theory can be used to calculate extinction cross section (aj, asymmetric parameter (g), 
and the single scattering albedo (o>) for a known refractive index (which depends on 
wavelength and composition) and particle radius; (2) the particle size distribution, n(r), 
is uniform within each pixel; and (3) the volcanic clouds are continuous, i.e. A =1. Based 
on assumption (3), (1) can be expressed as: 

I i ~il-R i (r m ,X c ))B(T 9 )*t i lr m .X c ) (B(T S )-B(T C )) (11) 

where r e is the effective radius of the spherical particle for given size distribution n(r), i.e. 


J:jcr 3 n(r) dr 

(12) 

Jti r 2 n(r) dr 



and x c is the optical depth defined by 

x e -L* fQ« rt Jcr a n(r)dr (13) 

where L is the geometric thickness of the clouds, Q cn is the efficiency extinction factor 
calculated by Mie theory, and c=K7 2 Q en is extinction cross section. 

Once the simultaneous pairs of radiance (I it /j) are calculated for varying r e and x c , 
the corresponding brightness temperature pairs (T if 7]) can be obtained by the rearranged 

version of Plank’s formula: 



( 14 ) 


1 . 43879 x 10 * 

Ti * X^lnT [3 .74l51xl0 8 Xr/rtJ 1 ] +1) 


where \ is wavelength in microns. 

Theoretical calculations based on the model are shown in Figure 1, where the soUd 
lines represent the effective radius of volcanic ash particles and dashed lines represent the 
10.8 pm (band4) optical depth of the cloud. The model reveals a nonlinear relationship 
of temperature difference with brightness temperature due to the nonlinear relation of 
optical depth with wavelengths. The theoretical calculation shows that the lowest 
temperature difference (LTD) of band 4 and band 5 linearly depends on the temperature 
difference of the cloud top and the surface (T, - Tc), i.e., when T, -T, increase, LTD will 
linearly increase. If the surface temperature T, is 273"K and the cloud top is 213"K (the 
conditions of our test data set), LTD cannot be lower than -3CfC, and the brightness 
temperature corresponding to LTD is about 22TK. It is found that the characteristics of 
negative temperature difference will disappear completely if the particle size is greater 

than 5 pm (see fig. 1). 

In figure 1, let S(r # ) be a variable that records the size of area embraced by the 
straight line which presents the appearance of 0 degree temperature difference and the 
curved line which represent the change of the temperature pairs (T 4 , 7> T s ) along with 
varying degree of transparency at different effective radii r. So S is a function of r, A 
plot of normalized S (maximum S is assigned to 1) versus the particle radius (O shows 
in Figure 2. S(rJ is a parameter that presents the retrieval sensitivity of the effective 
radius. If S is a monotonic function, then the observed temperature pairs (T 4 , T 4 - T s ) can 
be uniquely related to effective radius and optical depth pairs (r c , O . In our example (fig. 
2), S(r e ) is a monotonically decreasing function in this interval from 0.8 to 4.3 pm. 


STUDY CASE 


As an initial test for our model we examined volcanic cloud data taken from the 
August 19, 1992 Crater Peak/Spun eruption . Nine digital images of AVHRR were 
received from the NOAA 11 and 12 polar orbiting satellites. Previous study (Schneider 
et al, 1993) indicates that the images taken early in the eruption are apparently very rich 
in water droplets and/or ice and of large particle size. This makes the spectral signal 
similar to a meteorological cloud. However, when the drifting volcanic clouds dry out and 
the particle siz e becomes small during transport and dispersion, the clouds take on spectral 
properties dominated by fine volcanic ash. For the purpose of this study, the selected data 
was taken at 1338 GMT on August 19, 1992, about 13 hours after the onset of the 
eruption and 9 hours after it ended. At this point the cloud was located over the Gulf of 
Alaska, more than 300 km from Spun - Volcano. We used a sample frame of 1.1 km 
resolution Local Area Coverage (LAC) data of cloud surface temperatures of band 4 and 
band 5. In the sample frame, composed of 150 x 100 pixels, covering an area of about 
18150 km 2 , the volcanic clouds overlapped low-level meteorological clouds (Figure 3). 

Parameter selection 

In this study, we don’t have refractive index directly measured from Spurr ash. We 
assume that the volcanic clouds contain only volcanic ash and we used refractive index 
data on ash obtained by Volz (1973) and Pollack, Toon et al.. (1973). The six samples 
provide a good variety of volcanic ash types with crystalline andesite, crystalline basalt, 
glassy basalt, and glassy rhyolites (table2). The composition of Crater Peak/Spurr eruption 
is andesite, similar to sample 1, therefore we focus the refractive index of sample 1 and 
test the sensitivity of refractive index by comparing results using the other samples. 

The cloud top temperature (TJ is chosen at 213°K, which is determined by the coldest 
part of the earlier images of the volcanic clouds, where the absolute temperature difference 
between band 4 and band 5 is less than 0.5°C. The surface temperature (7\) is chosen 
273° K, which is determined from the band 4 and band 5 calibrations based on 
observations of areas, that surround the volcanic cloud in the AVHRR images, which are 



free of meteorological clouds, except 
underly the volcanic cloud. 


for the homogeneous low deck of clouds that 


Retrieval of particle sizes and optical depths 

Figure 4 shows results obtained from actual pixels in the sample frame superimposed 
on the calculated curves of figure 1. Most of points in the frame cluster between particle 
sizes of 1 and 4 pm (figure 4). To understand the spatial distribution of effecnve 
radiance, twenty theoretical curves from 0.5 to 5 pm with increments of 0.25 pm were 
calculated to analyze the AVHRR data in more detail. The minimum effective radius is 
to be 0.5 pm rather than 0 because of the existence of multiple solutions of our 
model for radius less than 0.8 pm (note the non-decreasing function S(rJ fraction of 
Figure 2). The frequency distribution shows pixels with an effective radius of less than 
0.8 pm are only 1.36% of the total pixels, and the mean radius and standard deviation are 
2.67 pm and 058, respectively (Figure 5). Mapping effective radii allows examination of 
the spatial distribution of particle size (Figure 6). The contour map shows that the 
effective particle sizes in the frame vary from 2.5 and 3 pm at the center, to between 3 
and 4 pm in the west, and to between 2 and 1.5 pm in the east There is a dramatic 
difference in particle size at the east and west edges of the frame (also the edges of die 
clouds), and the center is relatively homogenous (also the center of clouds). This variation 
could be due to changes in cloud altitude with dispersal as volcanic ash particles are being 
fractionated. A possible explanation consistent with the observed data is that the eastern 
edge of the cloud is higher and has moved faster than the western edge. It is also possible 
that at the edges of the cloud, the assumption of continuous coverage is violated, which 

gives us spurious results. 

The spatial distribution of optical depths of band 4 is shown in Figure 7, where twenty 
curves of optical depth with increments of 0.1 are calculated for retrieval. The optical 
depth in the sample test frame varies from 0.1 to 1.25. the mean value is 0.65 and 
standard deviation 0.28. The largest values of optical depth appear at the center of the 


clouds, and the smallest ones are at the cloud edges. 


Estimation of the total mass 

Using a density of 2.6g cm* 3 (Neal et al, in press) for the volcanic particles, the total 
mass of the frame is the accumulation of pixel-scale mass, Le. 

woo yiso L u.n s . p C*nr 3 n(r) dr 

a tOt Ml Z^i-1 1 r J j 




(i# ,/) ^ (i, j) 


J%r 2 n (r) dr 
jQ~iir^n{r) dr 


( ton) 


(15) 


where p(g/cm 3 ) is particle density, r e (ij) (pm), and t?* are retrieval effective radius and 
optical depth of pixel (ij), respectively. is geometric thickness of pixel (ij), and S is 
the area of each pixel. If the particle distribution n(r) reduces to uniform distribution, i.e. 
n(r)-l , then the inner term of the above sum simply becomes 1.21x4/3xpr/ v V ,J tyQ <I /r ( , M ). 
The £2e» at different effective radius are listed in table 3. 

Using the retrieval of r e (iJ) , ^ and Q ca discussed above, the total mass of volcanic 
ash in the study frame is about 36,200 tons. By extending the result to the whole cloud 
(about 100,000 km 2 in total area, see figure 4), the total mass is about 0.2xl0 6 tons, which 
is about 0.56% of the total volcanic ash measured in the deposited ash blanket (36X10 6 

tons) (Neal et al, in press). 

Since we have little direct information for the size distribution of the study cloud, the 
above estimation is simply based on the assumption of uniform size distribution. Our later 
work shows that there is no big change for the total mass due to several possible size 
distributions. This will be discussed in the next section. 


CONCLUSIONS AND DISCUSSIONS 

Theoretical study and application to a volcanic cloud imaged by AVHRR radiative 
transfer calculations in the infrared window channels can be used to retrieve the effective 


radius of volcanic ash particles and optical depths of clouds from AVHRR multispectral 

images. The major conclusions from this study are: 

(1) Volcanic clouds have negative brightness temperature differences (band 4 - band 
5) only if they have a dominance of particles with radius less than 5 pm. 

(2) Our model works best when there is large difference between the temperature of 
the underlying surface and the volcanic cloud. The lowest temperature difference (band 
4 - band 5) of the volcanic cloud is a linear function of the temperature difference 

between the underlying surface and the volcanic clouds (T, - TJ. 

(3) This method can be used to interpret volcanic clouds with a dominant effective 

radius between 0.8 and 4.3 pm. 

(4) The mean radius and the optical depths within the test frame of a 13 hour old 
August 1992 Crater Peak/spurr volcanic clouds are determined to be 2.67 pm and 0.65, 
respectively, based on the two-band model. The estimate of the mass of the volcanic cloud 
particles is about 28 - 36 x 10 3 tons in the frame, and about 0.15 - 0.2xl0 6 tons in the 
whole cloud, which is about 0.4 - 0.56% of the total volcanic ash measured in the ash 

blanket (36xl0 6 tons). 

The temperature difference model for a homogenous parallel cloud is completely 
controlled by the refractive index of particles in the cloud (which is dependent on their 
chemical composition, environmental variables T, and T c , the shape of the particles, and 
the size distribution of particles). Uncertainties in any of these parameters will influence 
the accuracy of the retrieval. At least for the Spurr cloud in our study, the model gives 
reasonable results (Figures 4-8). Additional studies of other volcanic clouds are needed 
to validate the model and to determine the conditions which favor accurate retrievals. 

Uncertainty about refractive index of volcanic ash 

In order to investigate the sensitivity of our model to uncertainty about refractive 
index, we compared results using the data from a variety of samples (figure 8; table 4). 


Some of the refractive index data (sample 4, 5 and 6) cause pixel points to fall outside of 
the fields of calculated results, we interpret this to reflect a lack of similarity in 
composition for those samples to the Spurr ash. For the three samples that do match spun: 
well, the effective radius determined ranges from 2 to 2.6 jim,the optical depth from 0.6 
to 0.67 and the mass in the frame from 28,000 to 36,000 tons. It would obviously be 
desirable to have refractive index data on the Spurr ash, but the sensitivity of this 
uncertainty is not too serious at this point. 

Silicates and sulfuric acid aerosols 

Silicates (volcanic ash) and sulfuric acid aerosol particles are both found in volcanic 
ash clouds (Rose et al., 1980; 1987; Tabazedeh and Turco, 1993). Prata (1989a) showed 
that either silicates or sulfuric acid can produce negative apparent temperature differences. 
We used refractive index data for sulfuric acid aerosols (Halperin and Murcray, 1987) to 
calculate model results analogous to those for silicates, for a variety of dilutions of H^SO*. 
Results show that aerosol with >50% K 1 S0 4 are indistinguishable from silicates with two- 
band AVHRR retrievals. The possible confusion of the two types of particles depends 
partly on the age of the volcanic cloud because ash is the dominant particle component 
of young ash clouds (Rose et al., 1980, 1987), but is removed from the atmosphere 
relatively quickly, in days to months (Turco et al., 1983; Pinto et al., 1989) while volcanic 
sulfuric acid may persist in the atmosphere for several years (Rampino et al., 1988; 
Bernard and Rose, 1989). For the Spurr cloud we used in our study, silicates probably 
exceed sulfuric acid particles by orders of magnitude. After atmospheric residence times 
of a week or more, confusion between the two particles types would be more of an issue. 
We have begun to investigate where three-band retrievals, involving band 3 of AVHRR 
(3.5 pm), can be used to distinguish between concentrated HjSO* and silicate components, 
it may also open the range of retrievable sizes, but the intensity of band 3 radiation in 
nighttime AVHRR may not be sufficient During the day it is useless due to reflected 

solar. 



Temperature difference between the volcanic cloud and the warmer subsurface 
The retrieval model we have developed is highly sensitive to the temperature 
difference (T, - TJ between the warmer surface (land, sea or clouds) and the top of the 
volcanic cloud. We expect drat (T. - Tj values will be typically positive, because 
tropospheric temperatures decrease rapidly with height and drifting volcanic cloud particles 
will equilibrate with the temperature of surrounding air. However the temperature of the 
surface can vary a lot, being highest on warm summer days in which there are no 
meteorological clouds (T. - T. up to as high as lOO’C) and lowest on days when high cold 
clouds underlie the volcanic cloud or in winter and night when die surface temperatures 
may be much lower. This effect is shown by comparing Figure 9 with Figure 1, which 
helps explain why some AVHRR images are such more successful in mapping and 
discriminating volcanic clouds (Schneider et al„ 1993) than others (Schneider and Rose, 

1993). 

Size distribution of particles 

The mass calculation is obviously dependent on the size distribution of particles (note 
the equations (13),(14),and (16)). The assumed uniform distribution is undoubtedly too 
simple to be reliable. Two other size distributions, Gamma and Lognormal, were 
suggested by many authors based on experiments and measurements (Deirmen, 1969; 
Prata, 1978). From the measurements of particles in the 1990 volcanic clouds of Mount 
Redoubt (Hobbs et al„ 1991), volcanic dust has Lognormal size distribution with 
parameters p=-1.36, 0=0.74 and mean radius=0.808 pm. We assumed o=4).74 in this 
study. Table 5 lists the comparison of the mass estimates based on different size 
distribution in the frame scale retrievaL It can be seen that the mass estimates vary only 
slightly from the uniform distribution for some possible gamma and lognormal 
distributions which could span the ranges of volcanic clouds, based on limited study of 
size distributions of particles by aircraft-based studies (Rose et aL 1980; Hobbs et al, 

1991). 

The table 5 implies the following conclusion; (1) when the variance increases in 


size distributions, the estimated mass will increase, as well ; (2) when lognormal 
distribution is chosen, an effective radius of 2 to 2.6 pm equals to a lognormal average 
of particle radius of between 0.7 and 0.9 pm, which coincidentally matches the average 
radius (0.808 pm) of airborne measurement of 1990 Mount Redoubt eruption cloud 
(Hobbs et aL, 1991) ; (3) when six refractive index samples are considered in lognormal 
size distribution, the estimated mass in the whole cloud is in the range of 280 - 290X10 3 , 
i.e. about 0.78-0.80% of the total volcanic ash measured in the ash blanket The smaller 
change of mass within the samples than within the distributions indicates that the size 
distribution which is used is more important than composition of the volcanic dust for 

mass estimation. 

Spatial resolution 

The method used in this study assumes that the volcanic cloud forms a well defined 
homogeneous single layer in each pixel. However, real volcanic clouds do not exactly 
meet this assumption. In fact, the particle size distribution in Figure 6 shows that cloud 
in the frame is not a single parallel layer, but probably has a range of altitudes. Therefore, 
the accuracy of retrieval will be improved by a higher spatial resolution of the remote 
sensor. A smaller sized pixel win reduce the effects of variable altitudes and the cloud 
within each pixel can be better approximated as a layer. Thus the use of LAC data (1.1 
km resolution) may be significantly better than GAC (4 km resolution) AVHRR data. 
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Table 2 . RefracUve index of different samples 


















Table 3. Efficiency extinction factors of 10.8 pm for volcanic ash 











Table 4. Pixel-scale retrival of masses for different samples 










Table 5. Frame scale retrieval for different size distributions 


Samples 


m (2) (3) (*) I & 

And^ite Basalt Glassy basalt Glassy rhyolite | Glassy Rhyolite 



mean effective radius | 2.6675 2.4736 


mean optical depth 11 0.6510 ^^6747^ 


2.0665 


a6057 


2.0252 


0.6399 


0.6379 


2.0027 


0.5929 


Uniform distribution: N(r^*l. meaner,,, vtr=0. 


mean pa rticle radius 
masa in a p ixel 
mass in the frame 
ma ss in the cloud 

mean particle radius 
mass in a pixel 
mass in the frame 
mass in the cloud 

mean particle radius 
mass i n a pixel 
mass in the frame 
mass in the cloud 


minimum radius 


maximum radius 


rnstimiM ndltK to get 
negative vtk» of T4-T5 


16675 

14009 

35.671 


14736 

11340 

31.705 


10665 

1.7851 

26,521 


10252 

1.9984 10027 

13978 

13677 16965 

35,624 

35,177 40,061 

196,000 

193,824 221,000 


I 197.000 1 175,000 146,000 I 196,000 | 193.824 j g 

Modified Gamma Dist ribution: n(rMN(rJ6) T /6!>*exp(-rj/6). mean=7r./9. var=7r. fSl 

|| 2.0747 1.9239 1.6073 | 1.3752 | 1S543 


1.9239 

1.6073 

1.5752 

1.5543 

14763 

11085 

2.6421 

16009 

36,790 

31,327 

39,253 

38,642 

203,000 

173,000 

216,000 

213,000 


2 fra? 1.9239 1-6073 1.5752 l£43 L5577_ 

2.6323 14763 11085 16421 16009 18647_ 

39,109 36,790 31,327 39,253 38 ' 642 42>56 1- 

215.000 203,000 173,000 216,000 213,000 | 235,000 

Lognormal Distribution: nfrW l/gwl^exptQnr-p)^}. var=e*" rt (e’ tf -l ). _ 

0,8922 | 0.8274 0.6912 ~~ 0.6774 0-6684 °- 6698 . 

3.4864 3.4387 18341 3.5151 3 ~ 4 ^ 3 ^ 3J582_ 

51,798 51,089 41106 51223 51 - 452 51864_ 

285.000 281,000 231000 388,000 283 ’ 000 291,00< ? 

The sensitive range of effective radius (pm) which can be uni quely retrial: 


0.8274 

0.6912 

0.6774 

0.6684 

0.6698 

3.4387 

18341 

3.5151 

3.4631 

3.5582 

51,089 

41106 

51223 

51.452 

51864 

281,000 

231000 

288,000 

283,000 

291,000 



Note: In the above listed number, the unit of mas, is in ton, and radius is micron, a in lognomul is chosen 0.74, and P-MO- 
2.5c 2 is the mean of mean of ln(r). 
























Figure Captions 


Fijuie 1 Two-band temperature difference model at 10.8 pm and 11.9 pm. The solid 
lines represent different effective radius, and the dashed lines represent the change 
optical depths at 10.8 pm along with particle radius. 

9 Effective radius and normalized negative temperature domain areas. The 
no ^ of temperature difference model to the effective 

radius. 

Figure 3. The AVHRR eruption cloud of August 19. 1992 Crater Pe^/Spurr at H38 
GMX The center square is the study frame, with an area of about 165 km x 1 10 . 

Eionre 4 Simulated temperature pairs, the temperature differences (10.8 pm - 12 pm) 
represent the observed data. 


Figure 5. Frequency distribution of effective radius. 


Fmure 6 Contour map, for the study area, of the effective radius in pixel-scale 
retrieval/lhe numbers in die contour map are the effective radius in microns. Numbers 
on the axes indicate distances among the edge of 1.1 km pixels. 


Figure 7. Contour map of the optical depth at 10.8 pm (band 4). 

Figure 8. The comparison of models for different refractive index. 

Figure 9. The two-band temperature difference model for T, - T=10°C. The absolute 
maximum temperature difference is reduced to about 3°C. 


Temperature difference (band4 - band5) ( C) 
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Temperature difference 




